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Employing the 12 Nal detectors in the Fermi GBM, the Earth Occultation Technique (EOT) can 
be used to measure the fluxes of x-ray and gamma-ray sources. Each time a source passes behind the 
Earth (or emerges from behind the Earth), a step-like feature is produced in the detector count rate. 
With a predefined catalog of source positions, the times of the occultation steps can be calculated, 
the individual steps fit, and the fluxes derived. However, in order to find new sources and generate a 
complete catalog, a method is needed for generating an image of the sky. An imaging algorithm has 
been developed to generate all-sky images using the GBM data. Here we present imaging results 
from ~ 2.5 years of data in the 12-25 keV and 100-300 keV energy bands. 



I. INTRODUCTION 

The hard x-ray/soft gamma-ray sky is extremely 
variable, and so the ability to continually monitor the 
sky at these energies is important. Especially above 
^200 keV, imaging with the wide-field coded aperture 
technique is difficult [1] so other detection techniques 
are required. With these two considerations in mind, 
the Earth Occultation Technique (EOT) applied 
12 Nal detectors from the Gamma-ray Burst Monitor 
(GBM) onboard Fermi is used to provide a the large 
field of view for continuous observations of the sky 
at energies 12-1000 keV. EOT has been used success- 
fully to study sources over this range of energies with 
the Burst and Transient Source Experiment (BATSE) 
on the Compton Gamma Ray Observatory (CGRO) 
0E] and with GBM 1 (i . Because EOT requires an 
input catalog of source positions to monitor sources, a 
new imaging approach is employed to discover sources 
not in the catalog. The imaging algorithm that has 
been developed is a tomographic method that uses 
projections of the sky over the course of the orbital 
precession period of the satellite ( ~ 53 days) to gen- 
erate an image. Results presented here include all-sky 
images for ~ 2.5 years of GBM data for the 12-25 keV 
and 100-300 keV broad CTIME energy bands and a 
table of sources that have been added to the EOT 
input catalog based on the imaging results. 

II. EARTH OCCULTATION TECHNIQUE 

As demonstrated with BATSE, monitoring known 
hard x-ray/soft gamma-ray sources is possible with 



simple non-imaging detectors with EOT. The count 
rate in a detector viewing the source decreases like a 
step function every time a source is occulted by the 
Earth. A similar feature, but with opposite polarity, 
occurs in the detector count rates when the source 
comes out of occultation. Occultation times can be 
predicted using the coordinates of the source and the 
spacecraft positions. The occultation time is defined 
to be the time at which the transmission through the 
atmosphere is 50% at 100 keV. In calculating a flux for 
the source of interest, a 4-minute window of 0.256-sec 
CTIME count-rate data, centered on the calculated 
occultation time, is fit to a quadratic background plus 
a source model for each source that occults during 
the window. The source model is constructed from 
combining an assumed source spectrum with an atmo- 
spheric transmission model convolved with the detec- 
tor response. For each energy channel, a scale factor 
is derived from the fit to the data and averaged over 
the detectors viewing the source. The weighted aver- 
age scale factor for the source of interest is multiplied 
by the model flux to determine the final flux [H [B] . 



III. EARTH OCCULTATION IMAGING 

Since EOT requires an input catalog with predeter- 
mined source positions, another method is needed to 
search for sources absent from the input catalog. An 
imaging method conceptually similar to x-ray com- 
puted tomography used in medical imaging has been 
developed. When a source occults, the Earth's limb 
can be projected onto the sky resulting in a "slice" 
on the sky over which counts are integrated. During 
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the orbital precession of the spacecraft (~ 53 days 
for Fermi), the angle between the source position and 
the plane of the satellite orbit changes. Thus over the 
course of an orbital precession period a range of pro- 
jections are sampled. The projections from rising and 
setting occultation steps can be combined to generate 
an image of the sky and localize a source. The angular 
resolution of this method is limited by the finite time 
required to go from full transmission to full attenua- 
tion for a source. The duration of a step for a source 
occulting in the plane of the satellite's orbit is ~ 8 
seconds. As the angle between the orbital plane and 
the source increases, the duration of the occultation 
step increases by ~ 1/ cos(/3), where j3 is the elevation 
angle between the source and the orbital plane of the 
satellite. The angular resolution is 360° xAt/P, where 
At is the occultation duration and P is the satellite's 
orbital period ( ~ 96 minutes) . The angular resolution 
ranges from - 0.5° at j3 = 0° to - 1.25° at = 66°. 
When j3 > 66°, occultation no longer occur. 

Instead of a catalog of known source positions, our 
imaging technique starts with a catalog of virtual 
sources covering te sky with spacing of ~ 0.25°, re- 
sulting in a catalog of ~ 660,000 positions on the sky. 
With these virtual source positions and the spacecraft 
position history, occultation times can be predicted. 
A 4-minute window of GBM CTIME counts data cen- 
tered on each occultation time is selected for the de- 
tectors that view the step at less that 75° from the 
detector normal. A first-order response correction is 
applied to the data from these detectors. To be able to 
sum rise and set steps together, occultation windows 
for rise steps are converted to set steps by rotating 
the data about the occultation time at the center of 
the window. A weighted average of all the windows 
for a virtual source is calculated for an orbital preces- 
sion period. The averaged window is rebinned from 
the standard CTIME resolution of 0.256 seconds to 
2.048-second resolution, then passed through a differ- 
ential filter of the form 

\^]=i+fa + fb _ ST^3=i-fa-fb 

u i — p J 

Jb 

where r« is the number of counts in bin j, f a is the 
inner boundary of the filter, and fb is the outer bound- 
ary [7]. The values of f a and fb used here are 3 bins 
and 8 bins, respectively. The differential filter gen- 
erates a dip at the occultation time while smoothing 
background features, where the amplitude of the dip 
is related to the intensity of the source. To calculate 
the amplitude, the window is fit in two sections. Bins 
within ±(2/ a + /&) of the occultation time are fit to a 
polynomial. The outer portions of the window are fit 
with a spline function and joined by a straight line. 
The amplitude of the virtual source is found by taking 
the difference between the two fits at the occultation 
time. 



Averaged Occultation Window for Sco X-1 for Aug 12, 200c 
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FIG. 1: Top: Averaged window for 1 day for Sco X-1. 
Bottom: The filtered data with the central portion fit to a 
polynomial (green), and the outer background portion fit 
to a spline function (red). 

The top panel in Figure 1 shows the averaged occul- 
tation window for one day for a point at the location 
of Sco X-1. The bottom panel shows the filtered win- 
dow with the polynomial fit in green and the spline 
fit in red. The dashed lines are the boundaries for the 
polynomial fit. 

A. Systematic Effects 

Four systematic effects are explicitly taken into ac- 
count in the processing: 1) high-spin rate data, 2) 
South Atlantic Analomy (SAA) passes, 3) source con- 
fusion with bright sources, and 4) virtual sources in 
which the rise and set steps of a bright source dom- 
inate the spline background fit. The first two effects 
are removed by making cuts on the data to remove 
times when those issues are present. High spin rate 
effects occur when the spacecraft slews, which leads 
to a rapidly varying background during the 4-minute 
window. The SAA pass effect comes from activation 
in the detectors during passes through the SAA and 
results in large spikes in the data, which again leads 
to a rapidly varying background. 

The last two systematic effects are caused by bright 
sources. The issue of source confusion arises from the 
ambiguity of where along the Earth's limb the mea- 
sured flux originates. Since the GBM detectors have 
no direct imaging capability, the flux for an occulta- 
tion time can be attributed to anywhere on the sky 
along the Earth's limb. Consequently, bright sources 
(e.g. Crab, Cyg X-1, Sco X-1) cause "X" patterns 
that can extend for tens of degrees. These "arms" 
make it difficult to determine the presence of fainter 
sources as they overwhelm any flux that comes from 
possible faint sources. To mitigate the effects of source 
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confusion, an algorithm has been developed to ignore 
occultation windows for virtual sources when bright 
sources occult close in time ( < 11 seconds) but far 
away on the sky. 

The final systematic effect involves virtual sources 
at roughly the same declination of a bright source and 
a few degrees away in right ascension. The results 
of the fits are large negative amplitudes for virtual 
sources. The virtual source is close enough on the 
sky to the bright source so that both rise and set oc- 
cultation steps for the bright source occur inside the 
4-minute window, but the two sources are far enough 
away such that when the window is filtered, the rise 
steps and the set steps of the bright source remain sep- 
arated. Consequently, the steps from the bright source 
are treated as the background level, which causes the 
calculated amplitude of the virtual source to be ex- 
tremely negative as the background level is larger then 
the source at the occultation time. To reduce this is- 
sue the fit boundaries are shifted towards the occulta- 
tion time at the center of the window and reduce the 
effect of the bright source on the background fit. 

Figure 2 demonstrates this problem by using a 
source close to the Crab (a = 83.63°, S = 22.02°) 
at a = 80.75°, S = 22.0° for a precession period. The 
top panel in Figure 2 shows the rise and set occulta- 
tion steps for the Crab, roughly at bin 35 and bin 75. 
The bottom panel shows how the filtered Crab steps 
misrepresent the background and thus the calculated 
amplitude of the virtual source. 



Averaged Occultation Window for Crab for Aug 12, 2008 to Oct 4, 2003 
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Filtered Averaged Window for Crab for Aug 12, 2008 to Oct 4, 2008 
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FIG. 2: Top: Averaged window for 1 day for the Crab. 
Bottom: The filtered data with the central portion fit to a 
polynomial (green), and the outer background portion fit 
to a spline function (red). 



IV. RESULTS 

This imaging algorithm has been applied to 17 pre- 
cession periods of data, ~ 2.5 years, to generate all- 

TABLE I: Sources found with Earth Occultation Imaging 



Source Name 


Flux (mCrab) 


Significance 


3A 1822-371 


38.7 


26.1 


3C 273 


16.0 


16.8 


LMC X-3 


16.5 


14.0 


NGC 3227 


11.1 


13.4 


NGC 5506 


12.0 


12.7 


IC 4329A 


9.8 


10.9 


IGR J21247+5058 


13.1 


10.7 


NGC 5252 


9.8 


10.7 


Coma Cluster 


9.6 


10.7 



sky images in the 12-25 keV and 100-300 keV energy 
bands. The data cover August 12, 2008 to Febru- 
ary 12, 2011. Figure 3 and Figure 4 show the 12-25 
keV and 100-300 keV images, respectively. These are 
significance maps with green contours at 3.3, 5, 10, 
and 20o\ The sources in red denote sources previ- 
ously in the GBM EOT input catalog, while those in 
blue denote sources added to the input catalog based 
on the imaging analysis. These sources were iden- 
tified through cross-correlating features in the maps 
with sources in the Swift /BAT, INTEGRAL/SP1, and 
Fermi /LAT catalogs, which were then analyzed with 
EOT. A source is added to the input catalog if the 12- 
50 keV significance from EOT is at least 10tr. Table 
1 lists the sources found, their flux in the 12-50 keV 
band, and their statistical significance. 
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FIG. 3: All-sky image in the 12-25 keV band for August 12, 2008 to February 12, 2011 
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FIG. 4: All-sky image in the 100-300 keV band for August 12, 2008 to February 12, 2011 
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